Quantum dynamics of an atom orbiting around an optical nanofiber 



(N 

o 

(N 
Oh 

in 



^— > ■ 

=5 



>: 

in, 

00' 

(n: 

OV 

o: 



X 



Fam Le Kien, 1,2, EI K. Hakuta, 1 P. Schneeweiss, 2 and A. Rauschenbeutel 2 

1 Center for Photonic Innovations, University of Electro-Communications, Chofu, Tokyo 182-8585, Japan 
2 Vienna Center for Quantum Science and Technology, 
TV Wien - Atominstitut, Stadionallee 2, 1020 Vienna, Austria 
(Dated: September 14, 2012) 

We propose a novel platform for the investigation of quantum wave packet dynamics, offering a 
complementary approach to existing theoretical models and experimental systems. It relies on laser- 
cooled neutral atoms which orbit around an optical nanofiber in an optical potential produced by 
a red-detuned guided light field. We show that the atomic center-of-mass motion exhibits genuine 
quantum effects like collapse and revival of the atomic wave packet. As distinctive advantages, 
our approach features a tunable dispersion relation as well as straightforward readout for the wave 
packet dynamics and can be implemented using existing quantum optics techniques. 
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Since the earliest days of quantum mechanics, the 
study of localized, time-dependent solutions to bound- 
state problems has attracted considerable attention 
Quantum mechanical objects show both particle-like and 
wave-like behavior and this duality can be readily ex- 
plored by analyzing the dynamics of wave packets. Such 
localized states behave like classical objects and follow 
classical trajectories as long as they do not disperse. 
However, for nonlinear dispersion relations, the wave 
packet will spread out and the dynamics cannot be de- 
scribed anymore using classical physics. One of the most 
prominent witnesses of genuine quantum mechanical dy- 
namics is when the collapse of the wave packet is fol- 
lowed by a revival, i.e., a relocalization of the wave func- 
tion [2|, |3| . This intriguing quantum effect has been ex- 
perimentally observed in a number of systems, see Q and 
references therein, and dispersion engineering of wave 
packets as well as their collapse and revival dynamics 
constitute an active field of current research [4J-|6|. 

Here, we propose a novel experimental platform for 
the investigation of quantum wave packet dynamics, of- 
fering a complementary approach to existing theoretical 
models and experimental systems. Our approach fea- 
tures a tunable dispersion relation as well as straightfor- 
ward readout for the wave packet dynamics and can be 
implemented using existing quantum optics techniques. 
It relies on laser-cooled neutral atoms which are inter- 
faced with an optical nanofiber as initially proposed in 
Ref. Q and thoroughly analyzed in Ref. Q. In this sce- 
nario, an atom orbiting around the nanofiber can show 
quantum mechanical dynamics including an initially clas- 
sical orbiting motion of the atomic wave packet around 
the nanofiber, a spread of the atomic wave function (col- 
lapse) , a partial relocalization of the atom (fractional re- 
vival) [9] , and, finally, a full quantum revival of the orig- 
inal wave packet. In the following, we study this dynam- 
ics quantitatively, give analytical expressions for relevant 
time scales, and derive the functional dependence of the 
expected signals when absorptively probing the dynamics 
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FIG. 1. Schematic of trapping an atom in an orbit around an 
optical nanofiber by the evanescent wave of a quasi-circularly 
polarized red-detuned guided light field. 



using a nanofiber-guided, resonant light field. 

Consider an atom moving in vacuum around a silica 
nanofiber in a cylindrically symmetric potential U(r), see 
Fig. [TJ We use cylindrical coordinates {r,ip,z}, with z 
being the axis of the nanofiber. Due to the cylindrical 
symmetry of the system, the angular momentum L z — 
—ihd/dip of the atom with respect to the fiber axis z is 
conserved. In the eigenstate problem, we have L z = hm, 
where % is the reduced Planck constant and m is an in- 
teger, called the azimuthal quantum number. Hence, the 
atom's center-of-mass (COM) motional eigenstates can 
be written as ^ vmK = (2ir)- 1 R vm (r)e imv e iKz . Here, 
v is the radial vibrational quantum number, K is the 
wave number of the matter wave along the z direction, 
and R vm (r) is the radial part of the wave function. We 
perform the transformation R vm {r) = u vm [r)/\/r. The 
function u vm {r) is determined by the equation 
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where M is the atomic mass, U^\r) = U^f(r) + U(r) 
is the effective potential, with the centrifugal potential 
U^\r) = n 2 (m 2 -l/4)/(2A/r 2 ), while £ um is the energy 
eigenvalue for the COM motion of the atom transverse 
to the fiber. Thus, the radial motion of the atom orbit- 
ing around the nanofiber can be reduced to the motion 
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of a particle in the one-dimensional effective potential 

There exist stable bound states for the atom if the ef- 
fective potential U^g has a local minimum outside of 
the fiber |8(. This may happen if U is attractive, oppo- 
site to the centrifugal potential . In order to produce 
a cylindrically symmetric attractive potential, we send a 
circularly polarized red-detuned optical field of frequency 
u> through the nanofiber. When the nanofiber radius a is 
small enough, the fiber can support only the fundamen- 
tal mode HEn which generates an evanescent- wave 
guided light field around the nanofiber. 

The optical potential for the atom is given by U op t = 
— a(w)|E| 2 /4, where a{uS) is the real part of the atomic 
polarizability at the optical frequency u and E is the 
positive-frequency component of the electric part of the 
guided light field. The details of the calculations of the 
optical potential U op t for a cesium atom in the vicinity 
of a nanofiber can be found in Refs. We also take 

into account the attractive fiber-induced van der Waals 
potential t/ v dw 01 the atom d, EH- The combination 
of the optical potential and the van der Waals potential 
yields the potential U = U opt + J7 v dW- 

We plot in Figs. [2ja)-(c) the effective potential U^g 
and the ground-vibrational-state eigenfunction u m = 
u v —x >m for a cesium atom in its electronic ground state or- 
biting around the nanofiber with the azimuthal quantum 
number m. In the following, the fiber radius is 200 nm 
and the wavelength and power of the trapping light are 
1064 nm and 20 mW, respectively. We note that these pa- 
rameters are well accessible experimentally 12|, [l3| . Our 
calculations show that a trapping potential with a mini- 
mum outside the fiber can be formed when the azimuthal 
quantum number m is in the range from 430 to 530. We 
find that, in this region, the ground-vibrational-state en- 
ergy £ m = £ v —i im increases with increasing m and ex- 
hibits a small negative curvature, see Fig. Et^d). We note 
that we can tune the dispersion relation by varying the 
power of the guided light field, its wavelength, and/or 
the fiber radius. 

In the following, the motion of the atom along the fiber 
axis is disregarded because it is independent of the mo- 
tion in the fiber transverse plane and, consequently, in- 
dependent of the azimuthal quantum number m. Let 
\tp m ) be the ground state of the atomic COM motion 
transverse to the fiber with the angular momentum Tim. 
The wave function of this angular momentum state is 
V>m(r) = (2-Kr)~ 1 l 2 u m (r)e tm ' p and the corresponding en- 
ergy is £ m . Here, we have introduced the notation 
r = {r, ip}. Consider a wave packet \tp) = I^ m c m |V'm), 
which is a linear superposition of the angular momen- 
tum states \ip m ) with the corresponding probability am- 
plitudes c m . The temporal evolution of this superposition 
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FIG. 2. Effective potential Uffi (solid blue lines) and ground- 
vibrational-state eigenfunction u m (dashed red lines) for a 
ground-state cesium atom orbiting with the azimuthal quan- 
tum number m = 446 (a), 468 (b), and 510 (c) around the 
nanofiber. (d) Dispersion relation of the ground-vibrational- 
state energy S m . 



state is given by the wave function 



V>(r,t)=^< 



-iS m t/h 



"0m (r). 
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We assume that the amplitudes of the individual angular 
momentum states are of the standard Gaussian distribu- 
tion form, c m = (27rAm 2 ) -1 / 4 exp[— (m — mo) 2 /4Am 2 ], 
with a peak at m-o and a standard deviation Am, and are 
truncated at m mm < mo and m max > tuq. In our numer- 
ical calculations, we use mo = 468, Am = 6, m m in = 446, 
and m max = 510. 

Since we have mo 3> Am 3> 1, we can expand £ m in a 
Taylor series around mo according to 



£ mo + £' mn (m - m ) + £'' (m - m ) 2 /2 , (3) 



where £' m = d£ m /dm and £ m = d £ m /dm . In our case, 
we have \£^ \ <C £' m which leads to the existence of two 
different time scales. Inserting Eq. ([3]) into Eq. @ then 
allows one to identify these relevant time scales of the 
evolution of the wave packet 0] . The classical period of 
the rotation of the wave packet in the fiber transverse 
plane is given by T ro t = 2Trh/£!^ l0 . The characteristic 
collapse time is defined as the time at which the spread of 
phase differences between the various oscillatory terms in 
Eq. |2]) is about n, that is, when the interference is most 
destructive. This characteristic time is given by T co n = 
2y^U/ (\£ mg | Am). The revival time for the atomic wave 
packet is defined as the time at which the phase difference 
between two neighboring terms in Eq. ^ is 2-7T, that is, 
when the interference is most constructive. It is given 
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FIG. 3. Atomic probability density \ip(r, t)\ 2 in the fiber 
transverse plane at different evolution times. \ip{r,t = 0)| 2 
is a Gaussian wave packet of angular momentum states with 
mo = 468, Am = 6, m m j n = 446, and m max = 510. 



by T r cv = 4cirK/\£m \. For the parameters used in our 
numerical calculations, we have £' mo /2Trh = 214 kHz and 
|££J/2tt7i 2.52 kHz, which lead to T rot = 4.67 lis, 
Tcoii = 37.3 /xs, and T rcv = 794 /is. 

We plot in Fig. [3] the spatial profiles of the atomic 
probability density \ip(r, t)\ 2 in the fiber transverse plane 
for different evolution times. Initially, at t = 0, the wave 
packet is well localized, see Fig. 02[a). We observe from 
Fig. [3jb) that for t ~ 20 /xs, the wave packet has already 
spread significantly and | -0 (r , | 2 is delocalized along a 
circle in the fiber transverse plane. Figures |3tc) and[3Jd) 
show that, when the evolution time t is about 100 lis or 
200 lis, the wave packet partially relocalizes in the form 
of 4 and 2 sub-packets, respectively. These regions are 
the regions of T rev /8 and T rov /4 fractional revivals 
According to Figs. EJe) and|3jf), the probability density 
|^(r, t)\ 2 reforms into a structure with a single dominant 
peak when the evolution time t is about 400 lis or 800 lis. 
The structure in Fig. [3fe) is the result of the half revival 
realized at T rcv /2. Near this time, the wave packet re- 
forms with the original periodicity, but the phase of its 
orbiting motion differs from the initial wave packet 0. 
The structure in Fig. |3]Jf) is the result of the full revival 
at T rov . Near this time, the wave packet reforms with 



the original periodicity. Its peak value is reduced and 
its spread is increased with respect to the original wave 
packet due to cubic and higher order terms in the disper- 
sion relation. In the special case where T iev /T I0 t is an 
integer, the full revival is in phase with the initial time 
development 

In order to experimentally reveal the predicted wave 
packet dynamics, we propose to probe the orbiting atom 
by a weak, resonant, quasi-linear ly polarized guided field 
E p (r) of frequency uj p . Neglecting the effect of the fiber 
on the spontaneous emission rate and the detuning [Tij . 
the rate of the scattering from the atom is proportional 
to the overlap between the atomic wave packet and the 
probe field, 



7scaW« / |^(r,t)| 2 |E p (r)| 2 dr 



(4) 



Assume that the main axis of the polarization of the 
probe guided field is aligned at the azimuthal angle ip = 
0. For the following, we take advantage of the fact that 
the intensity of the quasi-linearly polarized fundamental- 
mode guided probe field exhibits an azimuthal modula- 
tion according to 15] 



|E p (r)| 2 ex 
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where , 



Zip , 



and ez v ' are the cylindrical compo- 



nents of the mode profile vector function [10|, Il5| . 

Hence, we find 
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\V m cos (2A£ m t/h) , (6) 



where A£ m 
and V m are 



(£m+i — £m-i)/2 and the coefficients B 
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)dr. (7) 



We plot in Fig. [4] the time dependence of -f sca {t)- Fig- 
ure |U[a) shows that 7 sca (i) oscillates with an initial vis- 
ibility of almost 40%. The oscillations result from the 
classical rotation of the atomic wave packet around the 
nanofiber. With time, the modulation amplitude reduces 
and 7sca(i) reaches a quasi-stationary value for t ~ 15 lis. 
The period of oscillations of 7 sca (i) as obtained from the 
numerical evaluation is X^fc^ — 2.33 /is, one half of the 
classical atomic rotation period T rot = 4.67 [is. This 
reduction is due to the fact that the intensity of the 
probe field is symmetric with respect to the reflection 
{r, p} {r, tp + 7r} in the fiber transverse plane, see 
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FIG. 4. Time dependence of the scattering rate 7 SC a(i) in the 
ranges from to 50 /is (a), from to 1000 ps (b), and from 
650 to 900 fis (c). The polarization of the guided probe field 
is aligned at the azimuthal angle ip = 0. 



Eq. ([5]). In Fig. [4jb) , clear-cut resumptions of the os- 
cillation of the scattering rate 7 sca (i) appear when the 
evolution time is about 200 [is, 400 /is, 600 /is, and 
800 fjs, corresponding to T rcv /4, T rcv /2, 3T rcv /4, and 
T rov , respectively. As a specific characteristic of our prob- 
ing scheme, the fractional revival of the wave packet at 
Trev/8, see Fig. EJc), and the fractional revivals at odd 
multiples of T rev /8 as well as all higher order fractional 
revivals Q do not give rise to a modulation of 7 SC a(*)- 
This can be readily understood considering the four-fold 
azimuthal symmetry of \ip[r, t)\ for these fractional re- 
vivals. In conjunction with the cos 2 (^-dependence of the 
probe field intensity, one can easily show that the mod- 
ulation amplitude of 7 sca (^) in Eq. Q vanishes. Figure 
Hfc) shows a zoom of the signal around T rcv . The visi- 
bility of the scattering signal is reduced to about 1/3 of 
the initial visibility, revealing the increased spread of the 
wave packet upon revival, see Fig. EJf). 

In order to further link the wave packet dynamics with 
the predicted signal, j S ca,{t), we expand A£ m in Eq. ([6]) 
into a Taylor series of the second order around the central 
azimuthal quantum number toq and find 



(8) 



The first term in Eq. (jHJ leads to the oscillations of 
the scattering rate with the period Tisc^ = nh/S^ — 
T ro t/2. This period is the same for the initial and re- 
sumed oscillations of the scattering rate. Note that 
this also holds for the fractional revivals at T rev /4 and 
3T rev /4 where two diametric wave packets orbit around 
the nanofiber. The second term in Eq. (j8j) leads to the 
falloff and the resumptions of oscillations in the scat- 
tering rate. We define the falloff time of r y sca (t) as the 
time at which the spread of phase differences between 



various oscillatory terms in Eq. (J6j) is about tt and find 

V^Tcou/A. For the resump- 
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L fall 



tion time of 7 SC a(*)j we obtain Tc^umc = 7T?j/|£^ | = 
Tev/4- Based on the parameters used in our numerical 



calculations, we find Tjfc — 2.33 /is, T^f* = 16.5 fis, 
and Tresumc — 199 /xs, in agreement with Fig.H) 

The observation of the predicted quantum dynamical 
effects for cold atoms orbiting in a light-induced poten- 
tial surrounding an optical nanofiber is within the scope 
of current nanofiber-based quantum optics experiments. 
In particular, loading of atoms and wave packet prepa- 
ration in stable orbits around the nanofiber should be 
possible by starting with a stationary trap in which the 
repulsion of the atoms from the nanofiber surface is ac- 
complished by means of a blue-detuned nanofiber-guided 
light field @, EM2I- 

After laser-cooling of the atoms 
into the vibrational ground state of the two-color opti- 
cal potential, the latter can then be set into rotation 
using polarization modulators while suitably turning off 
the blue-detuned field, eventually preparing the trap and 
wave packet proposed here. This system can be seen as 
a two-dimensional "artificial atom" where the orbiting 
atom and the nanofiber take the role of the valence elec- 
tron and the ion core, respectively, while the Coulomb at- 
traction is replaced by a tunable light-induced potential. 
In view of the tunability of the dispersion relation, this 
would then implement a versatile experimental platform 
for the investigation of quantum wave packet dynamics. 
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